The gene regulatory basis of bystander activation in CD8+ T cells.

CD8+ T cells are classically recognized as adaptive lymphocytes based on their ability to recognize specific foreign antigens and mount memory responses. However, recent studies indicate that some antigen-inexperienced CD8+ T cells can respond to innate cytokines alone in the absence of cognate T cell receptor stimulation, a phenomenon referred to as bystander activation. Here, we demonstrate that neonatal CD8+ T cells undergo a robust and diverse program of bystander activation, which corresponds to enhanced innate-like protection against unrelated pathogens. Using a multi-omics approach, we found that the ability of neonatal CD8+ T cells to respond to innate cytokines derives from their capacity to undergo rapid chromatin remodeling, resulting in the usage of a distinct set of enhancers and transcription factors typically found in innate-like T cells. We observed that the switch between innate and adaptive functions in the CD8+ T cell compartment is mediated by changes in the abundance of distinct subsets of cells. The innate CD8+ T cell subset that predominates in early life was also present in adult mice and humans. Our findings provide support for the layered immune hypothesis and indicate that the CD8+ T cell compartment is more functionally diverse than previously thought.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we perform single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we performed single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Integrated longitudinal multi-omics study identifies immune programs associated with COVID-19 severity and mortality in 1152 hospitalized participants.

Hospitalized COVID-19 patients exhibit diverse clinical outcomes, with some individuals diverging over time even though their initial disease severity appears similar. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity. In this study, we carried out deep immunophenotyping and conducted longitudinal multi-omics modeling integrating ten distinct assays on a total of 1,152 IMPACC participants and identified several immune cascades that were significant drivers of differential clinical outcomes. Increasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, NETosis, and T-cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma immunoglobulins and B cells, as well as dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to the failure of viral clearance in patients with fatal illness. Our longitudinal multi-omics profiling study revealed novel temporal coordination across diverse omics that potentially explain disease progression, providing insights that inform the targeted development of therapies for hospitalized COVID-19 patients, especially those critically ill.

Cyclone: an accessible pipeline to analyze, evaluate, and optimize multiparametric cytometry data.

In the past decade, high-dimensional single-cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation, which are computationally intense and difficult to evaluate and optimize. Here, we present Cytometry Clustering Optimization and Evaluation (Cyclone), an analysis pipeline integrating dimensionality reduction, clustering, evaluation, and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full-spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification but also enables the unsupervised identification of lymphocytes and mononuclear phagocyte subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on a variety of cytometry datasets, which will further power immunology research and provide a scaffold for biological discovery.

Distinct pulmonary and systemic effects of dexamethasone in severe COVID-19.

Dexamethasone is the standard of care for critically ill patients with COVID-19, but the mechanisms by which it decreases mortality and its immunological effects in this setting are not understood. We performed bulk and single-cell RNA sequencing of the lower respiratory tract and blood, and plasma cytokine profiling to study the effect of dexamethasone on systemic and pulmonary immune cells. We find decreased signatures of antigen presentation, T cell recruitment, and viral injury in patients treated with dexamethasone. We identify compartment- and cell- specific differences in the effect of dexamethasone in patients with severe COVID-19 that are reproducible in publicly available datasets. Our results highlight the importance of studying compartmentalized inflammation in critically ill patients.

Cyclone: an accessible pipeline to analyze, evaluate and optimize multiparametric cytometry data.

In the past decade, high-dimensional single cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation which are computationally intense and difficult to evaluate and optimize. Here we present Cyclone, an analysis pipeline integrating dimensionality reduction, clustering, evaluation and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification, but also enables the unsupervised identification of lymphocytes and mononuclear phagocytes subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on variety of cytometry datasets which will further power immunology research and provide a scaffold for biological discovery.

Astroglial toxicity promotes synaptic degeneration in the thalamocortical circuit in frontotemporal dementia with GRN mutations.

Mutations in the human progranulin (GRN) gene are a leading cause of frontotemporal lobar degeneration (FTLD). While previous studies implicate aberrant microglial activation as a disease-driving factor in neurodegeneration in the thalamocortical circuit in Grn-/- mice, the exact mechanism for neurodegeneration in FTLD-GRN remains unclear. By performing comparative single-cell transcriptomics in the thalamus and frontal cortex of Grn-/- mice and patients with FTLD-GRN, we have uncovered a highly conserved astroglial pathology characterized by upregulation of gap junction protein GJA1, water channel AQP4, and lipid-binding protein APOE, and downregulation of glutamate transporter SLC1A2 that promoted profound synaptic degeneration across the two species. This astroglial toxicity could be recapitulated in mouse astrocyte-neuron cocultures and by transplanting induced pluripotent stem cell-derived astrocytes to cortical organoids, where progranulin-deficient astrocytes promoted synaptic degeneration, neuronal stress, and TDP-43 proteinopathy. Together, these results reveal a previously unappreciated astroglial pathology as a potential key mechanism in neurodegeneration in FTLD-GRN.

Inhibition of MET Signaling with Ficlatuzumab in Combination with Chemotherapy in Refractory AML: Clinical Outcomes and High-Dimensional Analysis.

Acute myeloid leukemia patients refractory to induction therapy or relapsed within one year have poor outcomes. Autocrine production of hepatocyte growth factor by myeloid blasts drives leukemogenesis in pre-clinical models. A phase Ib trial evaluated ficlatuzumab, a first-in-class anti-HGF antibody, in combination with cytarabine in this high-risk population. Dose-limiting toxicities were not observed, and 20 mg/kg was established as the recommended phase II dose. The most frequent treatment-related adverse event was febrile neutropenia. Among 17 evaluable patients, the overall response rate was 53%, all complete remissions. Phospho-proteomic mass cytometry showed potent on-target suppression of p-MET after ficlatuzumab treatment and that attenuation of p-S6 was associated with clinical response. Multiplexed single cell RNA sequencing using prospectively acquired patient specimens identified interferon response genes as adverse predictive factors. The ficlatuzumab and cytarabine combination is well-tolerated with favorable efficacy. High-dimensional analyses at single-cell resolution represent promising approaches for identifying biomarkers of response and mechanisms of resistance in prospective clinical studies.

Robust partitioning of microRNA targets from downstream regulatory changes.

The biological impact of microRNAs (miRNAs) is determined by their targets, and robustly identifying direct miRNA targets remains challenging. Existing methods suffer from high false-positive rates and are unable to effectively differentiate direct miRNA targets from downstream regulatory changes. Here, we present an experimental and computational framework to deconvolute post-transcriptional and transcriptional changes using a combination of RNA-seq and PRO-seq. This novel approach allows us to systematically profile the regulatory impact of a miRNA. We refer to this approach as CARP: Combined Analysis of RNA-seq and PRO-seq. We apply CARP to multiple miRNAs and show that it robustly distinguishes direct targets from downstream changes, while greatly reducing false positives. We validate our approach using Argonaute eCLIP-seq and ribosome profiling, demonstrating that CARP defines a comprehensive repertoire of targets. Using this approach, we identify miRNA-specific activity of target sites within the open reading frame. Additionally, we show that CARP facilitates the dissection of complex changes in gene regulatory networks triggered by miRNAs and identification of transcription factors that mediate downstream regulatory changes. Given the robustness of the approach, CARP would be particularly suitable for dissecting miRNA regulatory networks in vivo.

In Vitro Evolution Reveals a Single Mutation as Sole Source of Src-Family Kinase C-Helix-out Inhibitor Resistance.

Understanding cancer cell drug resistance to protein-tyrosine kinase inhibitors, which often arises from acquired mutations in the target kinase, is central to the development of more durable therapies. Experimental systems that reveal potential paths to resistance for a given inhibitor and kinase target have an important role in preclinical development of kinase inhibitor drugs. Here, we employed a codon mutagenesis strategy to define the mutational landscape of acquired resistance in HCK, a member of the SRC tyrosine kinase family and therapeutic target in acute myeloid leukemia (AML). Using PCR-based saturation mutagenesis, we created a cDNA library designed to replace each codon in the HCK open reading frame with all possible codons. This HCK mutant library was used to transform Rat-2 fibroblasts, followed by selection for resistant colonies with A-419259, a pyrrolopyrimidine HCK inhibitor and drug lead for AML. X-ray crystallography has shown that A-419259 binding induces outward rotation of the kinase domain αC-helix, a conformation incompatible with phosphotransfer. Remarkably, only a single resistance mutation evolved during A-419259 selection: histidine substitution for threonine at the gatekeeper position in the kinase domain. Deep sequencing confirmed representation of nearly all other missense mutations across the entire HCK open reading frame. This observation suggests that A-419259 and other C-helix-out Src-family kinase inhibitors may have a narrow path to acquired resistance in the context of AML cases where Hck is an oncogenic driver.

Expression of myeloid Src-family kinases is associated with poor prognosis in AML and influences Flt3-ITD kinase inhibitor acquired resistance.

Unregulated protein-tyrosine kinase signaling is a common feature of AML, often involving mutations in Flt3 and overexpression of myeloid Src-family kinases (Hck, Fgr, Lyn). Here we show that high-level expression of these Src kinases predicts poor survival in a large cohort of AML patients. To test the therapeutic benefit of Flt3 and Src-family kinase inhibition, we used the pyrrolopyrimidine kinase inhibitor A-419259. This compound potently inhibits Hck, Fgr, and Lyn as well as Flt3 bearing an activating internal tandem duplication (ITD). Flt3-ITD expression sensitized human TF-1 myeloid cells to growth arrest by A-419259, supporting direct action on the Flt3-ITD kinase domain. Cells transformed with the Flt3-ITD mutants D835Y and F691L were resistant to A-419259, while co-expression of Hck or Fgr restored inhibitor sensitivity to Flt3-ITD D835Y. Conversely, Hck and Fgr mutants with engineered A-419259 resistance mutations decreased sensitivity of TF-1/Flt3-ITD cells. To investigate de novo resistance mechanisms, A-419259-resistant Flt3-ITD+ AML cell populations were derived via long-term dose escalation. Whole exome sequencing identified a distinct Flt3-ITD kinase domain mutation (N676S/T) among all A-419259 target kinases in each of six independent resistant cell populations. These studies show that Hck and Fgr expression influences inhibitor sensitivity and the pathway to acquired resistance in Flt3-ITD+ AML.

CDK2 kinase activity is a regulator of male germ cell fate.

The ability of men to remain fertile throughout their lives depends upon establishment of a spermatogonial stem cell (SSC) pool from gonocyte progenitors, and thereafter balancing SSC renewal versus terminal differentiation. Here, we report that precise regulation of the cell cycle is crucial for this balance. Whereas cyclin-dependent kinase 2 (Cdk2) is not necessary for mouse viability or gametogenesis stages prior to meiotic prophase I, mice bearing a deregulated allele (Cdk2Y15S ) are severely deficient in spermatogonial differentiation. This allele disrupts an inhibitory phosphorylation site (Tyr15) for the kinase WEE1. Remarkably, Cdk2Y15S/Y15S mice possess abnormal clusters of mitotically active SSC-like cells, but these are eventually removed by apoptosis after failing to differentiate properly. Analyses of lineage markers, germ cell proliferation over time, and single cell RNA-seq data revealed delayed and defective differentiation of gonocytes into SSCs. Biochemical and genetic data demonstrated that Cdk2Y15S is a gain-of-function allele causing elevated kinase activity, which underlies these differentiation defects. Our results demonstrate that precise regulation of CDK2 kinase activity in male germ cell development is crucial for the gonocyte-to-spermatogonia transition and long-term spermatogenic homeostasis.

The Notch signaling pathway promotes basophil responses during helminth-induced type 2 inflammation.

Type 2 inflammation drives the clearance of gastrointestinal helminth parasites, which infect over two billion people worldwide. Basophils are innate immune cells that support host-protective type 2 inflammation during murine infection with the helminth Trichuris muris However, the mechanisms required for basophil function and gene expression regulation in this context remain unclear. We show that during T. muris infection, basophils localized to the intestine and up-regulated Notch receptor expression, rendering them sensitive to Notch signals that rapidly regulate gene expression programs. In vitro, Notch inhibition limited basophil cytokine production in response to cytokine stimulation. Basophil-intrinsic Notch signaling was required for T. muris-elicited changes in genome-wide basophil transcriptional programs. Mice lacking basophil-intrinsic functional Notch signaling had impaired worm clearance, decreased intestinal type 2 inflammation, altered basophil localization in the intestine, and decreased CD4+ T helper 2 cell responses following infection. These findings demonstrate that Notch is required for basophil gene expression and effector function associated with helminth expulsion during type 2 inflammation.

MiR-146a wild-type 3' sequence identity is dispensable for proper innate immune function in vivo.

The prevailing model of microRNA function is that the "seed region" (nt 2-8) is sufficient to mediate target recognition and repression. However, numerous recent studies have challenged this model, either by demonstrating extensive 3' pairing between physically defined miRNA-mRNA pairs or by showing in Caenorhabditis elegans that disrupted 3' pairing can result in impaired function in vivo. To test the importance of miRNA 3' pairing in a mammalian system in vivo, we engineered a mutant murine mir-146a allele in which the 5' half of the mature microRNA retains its wild-type sequence, but the 3' half's sequence has been altered to robustly disrupt predicted pairing to this latter region. Mice homozygous or hemizygous for this mutant allele are phenotypically indistinguishable from wild-type controls and do not recapitulate any of the immunopathology previously described for mir-146a-null mice. Our results indicate that 3' pairing is dispensable for the established myeloid function of this key mammalian microRNA.

The Src family kinase Fgr is a transforming oncoprotein that functions independently of SH3-SH2 domain regulation.

Fgr is a member of the Src family of nonreceptor tyrosine kinases, which are overexpressed and constitutively active in many human cancers. Fgr expression is restricted to myeloid hematopoietic cells and is markedly increased in a subset of bone marrow samples from patients with acute myeloid leukemia (AML). Here, we investigated the oncogenic potential of Fgr using Rat-2 fibroblasts that do not express the kinase. Expression of either wild-type or regulatory tail-mutant constructs of Fgr promoted cellular transformation (inferred from colony formation in soft agar), which was accompanied by phosphorylation of the Fgr activation loop, suggesting that the kinase domain of Fgr functions independently of regulation by its noncatalytic SH3-SH2 region. Unlike other family members, recombinant Fgr was not activated by SH3-SH2 domain ligands. However, hydrogen-deuterium exchange mass spectrometry data suggested that the regulatory SH3 and SH2 domains packed against the back of the kinase domain in a Src-like manner. Sequence alignment showed that the activation loop of Fgr was distinct from that of all other Src family members, with proline rather than alanine at the +2 position relative to the activation loop tyrosine. Substitution of the activation loop of Fgr with the sequence from Src partially inhibited kinase activity and suppressed colony formation. Last, Fgr expression enhanced the sensitivity of human myeloid progenitor cells to the cytokine GM-CSF. Because its kinase domain is not sensitive to SH3-SH2-mediated control, simple overexpression of Fgr without mutation may contribute to oncogenic transformation in AML and other blood cancers.

Developmental Origin Governs CD8+ T Cell Fate Decisions during Infection.

Heterogeneity is a hallmark feature of the adaptive immune system in vertebrates. Following infection, naive T cells differentiate into various subsets of effector and memory T cells, which help to eliminate pathogens and maintain long-term immunity. The current model suggests there is a single lineage of naive T cells that give rise to different populations of effector and memory T cells depending on the type and amounts of stimulation they encounter during infection. Here, we have discovered that multiple sub-populations of cells exist in the naive CD8+ T cell pool that are distinguished by their developmental origin, unique transcriptional profiles, distinct chromatin landscapes, and different kinetics and phenotypes after microbial challenge. These data demonstrate that the naive CD8+ T cell pool is not as homogeneous as previously thought and offers a new framework for explaining the remarkable heterogeneity in the effector and memory T cell subsets that arise after infection.

Selective Inhibition of the Myeloid Src-Family Kinase Fgr Potently Suppresses AML Cell Growth in Vitro and in Vivo.

Acute myelogenous leukemia (AML) is the most common hematologic malignancy in adults and is often associated with constitutive tyrosine kinase signaling. These pathways involve the nonreceptor tyrosine kinases Fes, Syk, and the three Src-family kinases expressed in myeloid cells (Fgr, Hck, and Lyn). In this study, we report remarkable anti-AML efficacy of an N-phenylbenzamide kinase inhibitor, TL02-59. This compound potently suppressed the proliferation of bone marrow samples from 20 of 26 AML patients, with a striking correlation between inhibitor sensitivity and expression levels of the myeloid Src family kinases Fgr, Hck, and Lyn. No correlation was observed with Flt3 expression or mutational status, with the four most sensitive patient samples being wild-type for Flt3. Kinome-wide target specificity profiling coupled with in vitro kinase assays demonstrated a narrow overall target specificity profile for TL02-59, with picomolar potency against the myeloid Src-family member Fgr. In a mouse xenograft model of AML, oral administration of TL02-59 for 3 weeks at 10 mg/kg completely eliminated leukemic cells from the spleen and peripheral blood while significantly reducing bone marrow engraftment. These results identify Fgr as a previously unrecognized kinase inhibitor target in AML and TL02-59 as a possible lead compound for clinical development in AML cases that overexpress this kinase independent of Flt3 mutations.

Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq).

We provide a protocol for precision nuclear run-on sequencing (PRO-seq) and its variant, PRO-cap, which map the location of active RNA polymerases (PRO-seq) or transcription start sites (TSSs) (PRO-cap) genome-wide at high resolution. The density of RNA polymerases at a particular genomic locus directly reflects the level of nascent transcription at that region. Nuclei are isolated from cells and, under nuclear run-on conditions, transcriptionally engaged RNA polymerases incorporate one or, at most, a few biotin-labeled nucleotide triphosphates (biotin-NTPs) into the 3' end of nascent RNA. The biotin-labeled nascent RNA is used to prepare sequencing libraries, which are sequenced from the 3' end to provide high-resolution positional information for the RNA polymerases. PRO-seq provides much higher sensitivity than ChIP-seq, and it generates a much larger fraction of usable sequence reads than ChIP-seq or NET-seq (native elongating transcript sequencing). Similarly to NET-seq, PRO-seq maps the RNA polymerase at up to base-pair resolution with strand specificity, but unlike NET-seq it does not require immunoprecipitation. With the protocol provided here, PRO-seq (or PRO-cap) libraries for high-throughput sequencing can be generated in 4-5 working days. The method has been applied to human, mouse, Drosophila melanogaster and Caenorhabditis elegans cells and, with slight modifications, to yeast.

CTDB: An Integrated Chickpea Transcriptome Database for Functional and Applied Genomics.

Chickpea is an important grain legume used as a rich source of protein in human diet. The narrow genetic diversity and limited availability of genomic resources are the major constraints in implementing breeding strategies and biotechnological interventions for genetic enhancement of chickpea. We developed an integrated Chickpea Transcriptome Database (CTDB), which provides the comprehensive web interface for visualization and easy retrieval of transcriptome data in chickpea. The database features many tools for similarity search, functional annotation (putative function, PFAM domain and gene ontology) search and comparative gene expression analysis. The current release of CTDB (v2.0) hosts transcriptome datasets with high quality functional annotation from cultivated (desi and kabuli types) and wild chickpea. A catalog of transcription factor families and their expression profiles in chickpea are available in the database. The gene expression data have been integrated to study the expression profiles of chickpea transcripts in major tissues/organs and various stages of flower development. The utilities, such as similarity search, ortholog identification and comparative gene expression have also been implemented in the database to facilitate comparative genomic studies among different legumes and Arabidopsis. Furthermore, the CTDB represents a resource for the discovery of functional molecular markers (microsatellites and single nucleotide polymorphisms) between different chickpea types. We anticipate that integrated information content of this database will accelerate the functional and applied genomic research for improvement of chickpea. The CTDB web service is freely available at http://nipgr.res.in/ctdb.html.